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A fully antisymmetrized microscopic model is developed for light two-neutron halo 
nuclei using a hyper-spherical basis to describe halo regions. The many-body wavefunction 
is optimized variationally. The model is applied to 6 He bound by semi realistic Minnesota 
nucleon-nucleon forces. The two-neutron separation energy and the radius of the halo 
are reproduced in agreement with experiment. Antisymmetrization effects between 4 He 
and halo neutrons are found to be crucial for binding of 6 He. We also properly extract 
two-neutron overlap functions and find that there is a significant increase of 30%- 70% in 
their normalization due to microscopic effects as compared to the results of three-body 
models. 

O. 

1. Introduction 

In some light nuclei, the proximity to particle emission thresholds allows loosely bound 
nucleons to tunnel out into the classically forbidden region and form what is typically 
referred to as a halo. One neutron halos, such as n Be and 19 C, and two-neutron halos, 
such as 6 He and n Li, are the most common, although there are nuclear states with 
the halo formed by protons or more than two neutrons. All these systems have very 
small separation energies and unusually large matter radii, when compared to neighboring 
nuclei, and many of their properties are determined by the long-range part of the many- 
body wavefunction. Halo systems are not specific to nuclear physics; a review covering a 
number of fields can be found in [ [1] . In this work, we focus on two- neutron halo nuclei 
in which long-range Coulomb effects are not present. 

Traditionally, the theory of halo nuclei has been dominated by simple few-body models. 
These models place their faith on the fact that the valence particles forming the halo are 
partially decoupled from the rest of the system, the core. In this frame of mind, one 
typically assumes that the core degrees of freedom are completely frozen, while meticu- 
lously treating the relative motion between the valence nucleons and the core. Then, the 
many-body problem containing a two-neutron halo reduces to a core+n+n three-body 
one. For an overview of three-body techniques applied to halo nuclei, see for example [ 
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The crucial assumption of few-body models, a macroscopic inert core, is a mixed bless- 
ing: on one hand, it allows one to focus on the long-range core-valence correlations, on the 
other hand, it is undoubtedly a crude simplification of the many-body problem. Stemming 
from the inert core employed, the main drawbacks of three-body models are: i) unknown 
properties of the core, i.e. the internal details of the core are left out completely and, 
only when needed, they are provided in an ad hoc manner, ii) the improper treatment of 
antisymmetrization [|3], iii) the use of effective core-nucleon interactions constructed case 
by case in an ambiguous manner, and iv) the fact that even when the two-body interac- 
tions are fitted to (some) properties of all two-body subsystems, the resulting three-body 
separation energy is (much) smaller than the experimental value [OH]- This pathological 
lack of binding is often cured with an empirical three-body force or the two-body forces 
are refitted to deliver the right three-body binding. Furthermore, there are indications 
that for reaction calculations three-body wavefunctions may require additional renormal- 
ization to account for microscopic effects missing in the inert-core approximation [|S]. 
Some improvement over the inert-core picture is provided by few-body models incorpo- 
rating collective excitations into the core [ [71 [8] ; they, however, still suffer from the 
above-mentioned drawbacks of few-body models, and in addition are limited to applica- 
tions where the core is a good rotor or vibrator. Given all these drawbacks, the predictive 
power of three-body models of two-neutron halo nuclei is rather limited. 

The above-mentioned drawbacks of few-body models are eliminated in microscopic 
(cluster) models of light nuclei. Over the last decade, there have been tremendous ad- 
vances in brute force ab-initio methods and due to increasing computational power, these 
models have improved their accuracy and predictability. Amongst these, the Green's 
function Monte Carlo (GFMC) [El [TO], the no-core shell model [ dH [TJl E3! , and molec- 
ular dynamics models [ [HI UHl [16] have been successfully applied to light s- and p-shell 
nuclei. Somewhere between few-body and ab-initio models are microscopic cluster models 
in which some degrees of freedom are frozen to reduce the complexity of the many-body 
problem. Of particular relevance to this work is the stochastic variational method (SVM) 
[H7J, Q2] which inspired several aspects of the model here presented. Applications of SVM 
and some other cluster models to two-neutron halo nuclei can be found for example in [ 
HH [201 EH [221 [23] . 

Despite the success of existing microscopic (cluster) models applied to non-halo nuclei 
and their ability to reproduce basic properties of some halo species, such as (three-body) 
binding energies and radii, questions may arise about their ability to capture the long- 
range halo correlations. These correlations are important because the halo nucleons spend 
a considerable amount of time in regions distant from the core and, generally, it is this 
part of configuration space that contributes the most to reaction observables. As is known 
from three-body models, the wavefunction describing a two-neutron halo nucleus ought to 
fall off exponentially. Existing microscopic structure models do not pay close attention (if 
any) to asymptotic regions, instead they use the binding energy to assess the convergence 
towards an eigenstate. The convergence of the energy, however, does not necessarily 
guarantee the convergence of the wavefunction in long-range regions. Moreover, most 
microscopic models exploit computationally tractable bases, most commonly Gaussians of 
one sort or the other. In principle, it should be possible to capture the slower exponential 
decay by using a large Gaussian basis, but as argued in [[TJ, quality precedes quantity 



when it comes to halo nuclei; that is the shape of the basis functions matters more than 
the size of the basis. 

In addition, existing microscopic structure models fail to provide input to reaction 
calculations of two-neutron halo nuclei. To feed reaction calculations, by themselves for- 
mulated in a few-body picture, one would have to extract information about halo particles 
from a full microscopic wavefunction which is a non-trivial task. Even though recently we 
have witnessed some progress in this direction for two-body non-halo projectiles [[2H ES], 
most microscopic structure theories are still far from providing such few-body-like infor- 
mation about three-body-like halo nuclei. It is for this reason that in reaction calculations 
the structure of halo nuclei is taken from few-body models despite all their drawbacks. 

It is obvious that both few-body and microscopic structure models have their appealing 
aspects as well as drawbacks. It is the aim of this work to combine the best of the two 
worlds: to develop a fully microscopic model for light two-neutron halo nuclei bound 
by nucleon-nucleon interactions that would describe simultaneously short- and especially 
long-distance regions. The novelty of our model is the integration of few-body and ab- 
initio methods. We use a basis that describes short-range correlations and at the same time 
preserves the correct three-body asymptotics. This is achieved at a high computational 
cost since matrix elements cannot be evaluated analytically. Our ultimate goal is to 
provide microscopic structure information to be used in reaction calculations involving 
two-neutron halo nuclei. 

In this paper our model is applied to the simplest two-neutron halo nucleus, 6 He. 
Preliminary results for this case have been published in [ [26] . Since neither 5 He thought 
of as 4 He+n nor the n+n subsystem are bound, 6 He is a Borromean system [12] . To provide 
input for reaction calculations involving this nucleus, two-neutron overlap functions are 
properly computed from our microscopic model and, to our best knowledge, they are for 
the first time expressed in hyper-spherical coordinates. By doing so, these functions are 
directly applicable to some reaction calculations and they can be compared directly to 
three-body wavef unctions. 

This paper is organized as follows. In Section [2] the theoretical formulation of our 
model is developed. Section [3] outlines some technical details and numerical procedures 
employed. Results for 6 He are discussed in Section [U A summary can be found in 
Section [5] along with an outlook of possible future developments. 

2. Theory 

The model here presented has as its primary goal the applicability to reactions involving 
light two-neutron halo nuclei and thus needs to capture few-body long-distance features of 
these systems. A two-neutron halo nucleus will be described by an antisymmetrized prod- 
uct of a microscopic core and a valence part consisting of two neutrons, or schematically 
\l/ = A core ~ (core x valence). The terms "core" and "valence" refer to distinct pieces of 
the wavefunction prior to the action of the core-valence antisymmetrizer j\ core - val upon 
which nucleons from the two parts of the wavefunction become indistinguishable. At large 
distances, our wavefunction decouples into the three-body-like form \1/ — > core x n x n 
of a desired shape most naturally expressed in Jacobi coordinates. It was therefore our 
choice from the early start to express the entire wavefunction in Jacobi coordinates. As 



an additional benefit of using these coordinates, translational invariance is guaranteed. 
All details of the model can be found in [127]- 

More precisely, the microscopic core with a fixed total angular momentum and parity 
JJ ore , and isospin T core with projection M Tcore is described by an antisymmetrized wave- 
function $ j£ ore T core m t corresponding to the core's ground state. If desired, excited 
states of the core can be included in the future. The valence terms are drawn from a suit- 
able hyper-spherical basis. Each valence basis function ip carries total angular momentum 
J£ al , and isospin T va i with projection M TvaV which coupled to the core's quantum numbers 
give J 77 , T, and Mt for the full system. For neutron-rich two-neutron halo nuclei, such as 
6 He and n Li, T va i = 1 and Mx val = — 1, and the core-valence isospin coupling is trivial. 
It is understood that parities n carry the same subscripts as their corresponding J and 
that 7r = ir core ir va i. Then, the total wavefunction is written as: 

*J*Af J TM T = Y^ C J?ore F V al J£ al A"" ' " ' ^ [ $ ^ cor£ T core M Tcor£ ® A^ ^T val J^ l-l] j„ Mj y ^ • 

(1) 

By acting on valence particles only, the operator A val = Xa(~ 1) P -^ antisymmetrizes the 
valence part, whereas the core-valence antisymmetrizer j[ core - val = y^ ~ (— 1) P P p er _ 
mutes valence particles with those inside the core. In these operators, P are permutation 
operators, p are permutation parities, and A = A core + 2 is the mass number of the halo 
nucleus. The meaning of T va i is explained in Section 12.21 

The wavefunction in Eq. (jTJ) is constructed in two steps. First, the wavefunction of the 
core as a free nucleus is built within SVM, as briefly described in Section 12.11 Unlike many 
microscopic cluster models employing a simple Os-harmonic oscillator approximation to 
4 He, we use the best possible wavefunction for the core obtainable within SVM. By doing 
so, we hope to attenuate the problem of underbinding relative to the three-body threshold [ 
Q32 |2H]. Once the core wavefunction is optimized, its parameters remain unchanged in the 
subsequent minimization procedure. This implies that, as in other microscopic cluster 
models, distortion of the core due to valence neutrons is not accounted for explicitly, 
although some distortion is delivered through the core-valence antisymmetrization. In 
the second step, the valence part is constructed by drawing terms ipr val j^ i_i from a 
hyper-spherical basis described in Section 12.21 Valence basis functions contain discrete as 
well as non-linear continuous parameters that are optimized variationally, along with the 
linear coefficients c in Eq. (JTJ) , by minimizing the expectation value of the Hamiltonian 
H: 

(&\H\iff} 

E = -^y ( 2 ) 

as outlined in Section [3j To feed reaction calculations involving nuclei of interest, two- 
neutron overlap functions are extracted from the resulting optimized wavefunctions, as 
presented in Section 12.41 

2.1. The core in SVM 

Our choice of the microscopic model for the core has been motivated by the following 
factors. The model should provide accurate structure for the core; it needs to be extend- 
able to cores heavier than 4 He; and it must handle central and non-central nucleon-nucleon 





Figure 1. Sets of relative Jacobi coordinates x for a four-particle (pi,...,p 4 ) system: 
K-like (left) and H-like (right). Each Jacobi coordinate connects centers of masses of 
subgroups of particles. 



forces. Unlike for the valence part, there is no need to impose special asymptotic require- 
ments on the core. Finally, as mentioned earlier in this section, the wavefunction of the 
core should be expressed in Jacobi coordinates. It is for these reasons that the SVM model 
seemed to be the most appropriate. In this section, basic ingredients of SVM applied to 
4 He are outlined; more details on the method can be found in [[T7J [TBI 127] . 

In SVM, the core wavefunction is written as a linear combination of basis functions 6: 
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Jcore Tcore M T( , or 



A co ' e 2_^ c Vcorc j£ ore Tcorc M Tcore 4>r core j? ore T core M Tcor 



(3) 



with the projection of J core suppressed. The operator ¥ \ core antisymmetrizes particles 
inside the core. Basis states <ft take the form of correlated Gaussians expressed in Jacobi 
coordinates. For a four-particle system, there exist two different — K- and H-like — sets of 
Jacobi coordinates x = {xi,X2,xs} as shown in Figured] To accelerate the optimization 
of the core's wavefunction, both Jacobi sets enter Eq. ([3]) as they invoke different inter- 
particle correlations. For 4 He, a basis term in either Jacobi basis is given by: 



-L core ^ core ^-core ^T C or 



(£, A) = exp ( --xAx 1 [O h i 2 i 3 L 12 L(x) <8> Xs 12 s 12S s] J?ore t Ti2 t 123 t c , 



(4) 



The function 9i x i % i 3 l 12 l{x) is taken as a vector-coupled product of solid harmonics [ 

OhhhL l2L (x) = [[y h (x!) ® y h (x 2 )] Ll2 ® yi 3 (x 3 )] L . (5) 

The spin Xs 12 s 123 s an d isospin T Tl2Tl23TcoreMT parts consist of successively coupled single- 
particle spins and isospins: 



A,Sl2iSl23'S 



'Ti 2 T\ 2 sT core MT cor 



Sl2Z 



[[Xpt®Xp*]s„®Xpz]„ ®X P . 



[[r pi ® r P2 ] Ti2 ® r m ] 



'T 12 3®V 



T core MT cor 



(6) 

(7) 



For He, the Gaussian part in Eq. (J3J) contains a (3 x 3)-dimensional positive-definite, 
symmetric matrix A, specific to each basis term. The quadratic form xAx involves scalar 



products of Jacobi vectors: 

3 

JL-t\-JL = y J^-^jJu^ ' Jb j . I O ) 

Due to the symmetry requirement, the matrix A has only 6 independent elements and 
they are considered non-linear continuous variational parameters. Note that, although the 
Gaussian in Eq. (j4j) as a whole is a spherically symmetric object, it still carries angular 
information due to cross terms x, • Xj when the matrix A is non-diagonal as considered 
here. In such a case, numbers h,h, h in Eq. (EJ) loose their meaning of orbital momentum 
quantum numbers and can be treated as discrete variational parameters, instead. 

The composite index T core comprises other (quantum) numbers, elements of the matrix 
A, and the Jacobi channel identifier K or H, i.e. T core ={h,h,h, -^12, L, Si 2 , S123, 
S,Ti2,Ti23, K/H, A}. The sum in Eq. fl3]) was left without a summation index because 
in SVM the wavefunction $ is constructed term by term by minimizing the expectation 
value of energy with all components of T core optimized stochastically. Linear expansion 
coefficients c in Eq. fl3]) are determined via energy matrix diagonalization. 

2.2. The valence part in the hyper-spherical formalism 

To inspire the form of valence functions in Eq. ([T]), let us first consider a few-body 
approach in which a two-neutron halo nucleus is treated as a three-body core+n 1 +n2 
system. From among all three-body models, we adopt the formalism of [ [30] and [[31] pg. 
278]. For all necessary details, see [l2"7]. 

For a three-body core+ni+n2 system, one can define two different — Y- and T-like — sets 
of Jacobi coordinates x = {^i,^} as shown in Figure |2j These coordinates are further 
transformed into hyper- spherical coordinates {p,9,x} where x = {xi,£i} are spherical 
angles corresponding to vectors x. The hyper-radius p is related to the overall size of the 
system, while the hyper-angle 9 contains information about relative magnitudes of x\ and 
X2- In three-body models, the use of hyper-spherical coordinates is motivated by the facts 
that: i) the three-body Schrodinger equation reduces to a one- dimensional hyper-radial 
differential equation, and ii) more importantly for us, in the absence of long-range forces 
the three-body wavefunction of Borromean systems decays asymptotically for p — > 00 as: 

p~ 5/2 exp(-/tp), k 2 = 2m \E 3body \ /h 2 , (9) 

where E^ ody is the binding energy relative to the three-body core+n+n threshold, and m 
is the mass of a nucleon. In hyper-spherical coordinates, the volume element becomes: 

dV = (/ii/i 2 )" 3/2 p 5 sin 2 9 cos 2 9 dp d9 d Xl dx 2 , (10) 

where /ij, i = 1, 2 are dimensionless reduced mass factors corresponding to vectors x%. 

Upon the hyper-spherical transformation, the spatial part of the three-body wavefunc- 
tion can be separated into its hyper-radial, hyper- angular, and spherical parts, each of 
which can be conveniently expanded. The spherical part is written as a coupled product 
of spherical harmonics Yi irril (xi) and YJ 2m2 (x2). The hyper-angular part is expanded 
over eigenstates (p^ 2 (9) of the grand-angular operator appearing in the hyper-radial 
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Figure 2. Sets of relative Jacobi coordinates for a three-body core+n 1 -|-n 2 system: Y-like 
(left) and T-like (right). Each Jacobi coordinate connects centers of masses of subgroups 
of objects. These coordinates are not to be confused with those for the core in Figure [TJ 



Schrodinger equation. These functions are explicitly defined in terms of Jacobi polyno- 

l +— l 4-1 

mials P njac 2 ' 2 of the order rij ac = 0, 1, 2, . . . with hyper-momentum K = li + l 2 + 2n,j ac : 

#»(0) = N^ h (p 1 p 2 f /k sm h tfcos' 2 eP l n ^J M+ ^{cos2e). (11) 

The normalization constant N£ 2 is chosen to make these functions orthonormal with 
respect to a weight from the hyper-spherical volume element in Eq. ( fTUl) : 

f ip l k\0)ip l kH0) (/ii/i2)" 3/2 sm 2 ^cos 2 0d^ = 6 k ,k>. (12) 

JO 

Having in mind the importance of the asymptotic behavior of the wavefunction in 
Eq. 09]), the three-body model in [ [30] employs a hyper-radial basis from [ [32] with the 
desired exponential trend built into it: 



Here, l£ are associated Laguerre polynomials of the order ni ag . Note that this basis is 
just a suitable mathematical basis whose elements cannot be interpreted as hyper-radial 
eigenfunctions of the physical three-body system. Basis functions 1Z nia are orthonormal 
with respect to the weight factor p 5 from the volume element in Eq. (TTOj) : 



Kn lag (P, Po) Kn' lag (P, Po) P 5 dp = $n lag ,n' la 



(14) 



The completeness of the hyper-radial basis for any po > allows one to use any value 
of po, and yet reconstruct the asymptotically correct form of Eq. OH]) determined by an a 
priori unknown three-body binding energy. The index rii ag is independent from quantum 
numbers attached to spherical and hyper-angular parts of the wavefunction. 

Based on these three-body arguments, we have chosen for each valence term in Eq. (JTJ 
the following form: 

^r val j Zal 1-1 = K niag (p, po) % val J- al i-i(0, x), T val = {ni ag , Po, Jvai}, (15) 



where ^ is a generalized hyper-harmonic function in LS-coupling: 



% val jZ al i-i(9,Z) = vV 2 (8) 



[Y h {xi) ®Y h {x 2 )} L ®\x ni ® X 



"2J5 



(16) 



vat 



Here, \n t and r Ui are spinors and isospinors of valence neutrons. The composite index 
7„ a ; comprises all other numbers as well as the Jacobi channel identifier Y or T, i.e. 
■j va i = {K, li, l 2 , L, S, Y/T}. By construction, hyper-harmonics W corresponding to a 
given Jacobi basis are orthonormal in all components of 7„ a /. In the T Jacobi basis, the 
Pauli principle between valence neutrons is satisfied by requiring l\ + S = even; in the 
Y Jacobi basis, the exclusion principle is satisfied upon the action of A val in Eq. (JTJ). 
Valence angular momenta h, l 2 , L, and S are not to be confused with those for the core 
in Section 12.11 

2.3. The Hamiltonian 

In this work, the nuclear Hamiltonian includes kinetic energies Tj of all nucleons and 
two-body nucleon-nucleon potentials V#: 

A A 

No correction is needed for the kinetic energy of the total center of mass as the wavefunc- 
tion in Eq. (JTJ is expressed in relative Jacobi coordinates. 

In microscopic calculations, the choice of the effective nucleon-nucleon interaction is 
of crucial importance unless realistic forces are used. If a model is to have anything 
to do with the real physical problem, one must make sure that the inter-nucleon force 
is appropriate for all subsystems appearing in the model. In 6 He, valence neutrons are 
mostly in spin-singlet configurations [El EH]- The spin-singlet di-neutron state is unbound; 
however, many effective nucleon-nucleon interactions, such as the Volkov force [ [33] , do 
not distinguish this state from a bound spin-triplet neutron-proton state, the deuteron. 

In this work we use the semi-realistic Minnesota interaction [ 34J. This force reproduces 
the most important low-energy nucleon-nucleon scattering data and therefore it does not 
bind the di-neutron. The force renormalizes effects of the tensor force into its central 
component and binds the deuteron by the right amount assuming a proton and a neutron 
in a relative s-wave. It also gives realistic results for the bulk properties of nuclei in the 
lowest s-shell. Besides the central component, the potential contains spin, isospin, and 
spin-isospin exchange terms. The potential contains the mixture parameter u which can 
be tuned slightly to adjust the interaction strength. When supplemented by a spin-orbit 
force [[35], the Minnesota interaction reproduces low-energy a-nucleon scattering data. 
In [ [31] , it is advised to employ a short-range spin-orbit force; on that merit we use the 
spin-orbit parameter set IV from Table 1 in [ [35] . 

In Section HI we mostly show results obtained for the Minnesota plus the spin-orbit 
interaction (MN-SO). For comparative purposes, some calculations were performed in the 
absence of the spin-orbit force (MN). To reproduce the two- neutron separation energy 
0.97 MeV of 6 He, the mixture parameter u is changed from its default value 1.0 to 1.015 
(MN-SO) and 1.15 (MN). By doing so, we expect to obtain a realistic description of the 



halo even though the absolute binding of 4 He and 6 He is not reproduced. Essentially, the 
interaction mixture parameter is the only free parameter in our model. The Coulomb 
interaction is neglected as it should not affect the long-range behaviour of the neutron 
halo. 

2.4. Two-nucleon overlap functions 

Once the full wavefunction in Eq. (p} is optimized as outlined in Section [3j one can use it 
to calculate various observables. Although most of the observables we calculate are stan- 
dard, two-nucleon overlap functions deserve special attention. For general considerations 
regarding these functions, see [ [3S1 sect. 16.4.2] and [137]. 

The overlap integral between a two-neutron halo nucleus (\l/) described microscopically 
by Eq. flTJ and its own core in the ground state ($) is defined as: 



X. 



Mj ,1-1 

' J val 




e M Jcore T c°re M Ta 



I* 



J^MjTMt 



)■ 



The binomial factor accounts for the number of combinations to pick two out of A nucle- 
ons. The integration is done over all degrees of freedom in the core, and so the overlap 
integral X depends only on the degrees of freedom of two valencezl neutrons remaining 
outside the core. For neutron-rich two-neutron halo nuclei, the integral has a well defined 
isospin and its projection, 1 and -1, respectively, but it does not have a good angular mo- 
mentum. The integral can, however, be expanded in a complete set of hyper-harmonics 
of angular momentum J^ al introduced in Eq. ([TBI) : 



^Mjval 1 - X ~~ Z_^ C J C or J e M Jcore J val Mj ml G lval J* al (p) %val -C, 1 - 1 > 



(19) 



Ival ■J,,,, 



where C are Clebsch-Gordan coefficients. The expansion is carried out in the T Jacobi 
basis where the hyper-harmonics *3/ satisfy the Pauli principle between valence neutrons 
by construction (see Section I2~2"j) . The hyper-radial part in Eq. f[T9l is not expanded in 



the Laguerre basis from Eq. (TT3l) because the basis functions 7Z nia do not have physical 
significance. Instead, the overlap functions O are computed directly from: 



O-y ,-r (pO 

Ival -J,,„l vr ' 




, T CO re M Tcore ® % val J- ai l-l] j, Mj 



TM 7 



S(P - P') 



p° 



* 



J^MjT M T 



(20) 



with the integration carried over degrees of freedom of all nucleons. A meaningful cal- 
culation of overlap functions O requires both wavefunctions $ and \1/ to be normalized. 
Using overlap functions, a three-body-like core+n+n component of the wavefunction ^ 
can be written as: 



^.overlap 

* J^MjTMt 



Z2 ^"l J val [^JSore Tcore M Tcore ® % val J^l-l] 
7uai J va l 



J*MjT M T 



(21) 



2 In the overlap integral, the two distinct neutrons outside the core are called "valence" They are not 
to be confused with neutrons in the valence part of the wavefunction in Eq. (fj} where valence and core 
particles become indistinguishable upon the full antisymmetrization. 
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in a form analogous to that of the three-body wavefunction [[30]. Note, however, that 
the core $ in Eq. (1211) is fully microscopic, whereas the three-body wavefunction would 
contain an inert macroscopic core. Generally, the expansion in Eq. (j21j) would be a part 
of the fractional-parentage expansion of *£?. 

Overlap functions O satisfy a three-body-like Schrodinger equation with a source term 
[ |37] . Therefore, at least in the asymptotic region, wavefunctions from three-body models 
and our microscopically founded overlap function \fr overla P should behave similarly. On 
this merit, we can compare our results with those from three-body models at the level of 
wavefunctions rather than integrated observables. An overlap term characterized by ^ va i 
is referred to as an overlap or a three-body channel. 

For practical considerations, it is convenient to introduce modified overlap functions: 

u{p) = p 5 ' 2 O(p) (22) 

with a regular behavior near the origin. Again, in the absence of long-range forces in 
Borromean-like nuclei, overlap functions decay as: 

O(p) ^^ p- 5/2 exp(-/tp), u(p) ^^ exp(-Kp) (23) 

just like three-body wavefunctions in Eq. fl9]). Formally, the decay constant k in Eq. (|23|) 
is defined as in Eq. (Q, the main difference is that for overlap functions the three-body 
binding energy needs to be computed microscopically. For the Borromean nucleus 6 He: 

E 3body = £( 6 He) - £( 4 He) (24) 

where E( 6 He) and £'( 4 He) are binding energies corresponding to \l/ and $, respectively. 

Given the orthonormality of hyper-harmonics in Section I2.2[ the norms — spectroscopic 
factors — of overlap channels are given by: 



S y ,j* = Ol „ (p)p 5 dp= / ul „ (p)dp. 

>J U J U 



(25) 



In three-body models, spectroscopic factors give the probability of finding the system in 
a given channel ^ va i. 

3. Numerical details 

In Figure |3l we show the convergence of the binding energy of 4 He with the number 
of Gaussians included in the SVM wavefunction described in Section [2TT1 In converged 
MN and MN-SO (see Section I2.3P states containing 20 and 75 basis states, 4 He is bound 
by -30.85 MeV and -30.93 MeV, respectively. In both cases, all channels with / < 2 in 
Eq. ()5]) were present in the model space. 

The expectation value of energy in Eq. ([2]) involves multidimensional integrals which 
must be evaluated efficiently to perform a meaningful variational calculation. For ex- 
ample in SVM employing the basis of correlated Gaussians, many matrix elements can 
be evaluated analytically in closed form, and consequently SVM can afford a random 
trial-and-error variational search. In our model, however, the core is combined with a 
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Figure 3. Convergence of the binding energy of 4 He: MN (left) and MN-SO (right). 



functionally very different valence part, and these two become completely entangled upon 
antisymmetrization. In addition, there does not seem to be an analytical way of comput- 
ing matrix elements involving different core- valence permuted pieces of the wavefunction 
in Eq. (0Q). For these reasons, we are left with numerical evaluation of all matrix elements. 
For 6 He, the integrals in Eq. (J2J) involve 6 x 3 = 18 spatial and 2 x 6 = 12 spin-isospin 
dimensions. 

We use techniques of variational Monte Carlo (VMC) [ O [38J [39] to perform the vari- 
ational search. The integration space is sampled statistically to find the regions most 
relevant for a given physical problem. In nuclear physics, VMC has been successfully 
applied to the problem of light nuclei [ |9] ; however, the work here presented is novel in 
the choice of the trial wavefunction and therefore has its own challenges. In this section, 
we briefly mention the most important aspects of the optimization procedure, and advise 
the expert reader to [ |27] for a complete description. 

The greatest challenge we faced was to develop a robust and yet efficient method to 
navigate the parameter space. This aspect falls beyond the framework of VMC. We found 
that optimization methods used in atomic and molecular physics [ E0j HI] did not work 
adequately for the problem at hand, most probably due to the specific nature of non- 
central, state-dependent nuclear interactions. In addition, due to the closeness to the 
three-body break-up threshold, an arbitrary starting wavefunction is inevitably three- 
body unbound. As a consequence, since the trial wavefunction drives the scanning of the 
integration space, all optimization methods tend to break the nucleus into the core and 
two individual neutrons as long as the system is three-body unbound. The easiest way to 
solve this pathological problem is to constrain the radius of the nucleus while performing 
the energy minimization above the three-body threshold. In our model, this is achieved 
most readily by taking the same non-linear parameters po i n Eq. ( TT5|) for all valence terms. 
Once the system becomes three-body bound, po, still being the same in all valence terms, 
is adjusted to minimize the binding energy. 

Different optimization procedures were used for the two interaction cases, MN and 
MN-SO, of 6 He. In the absence of the spin-orbit force, the MN wavefunction contains 
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only spin-singlet valence terms. It is then possible to simply add hyper-angular and hyper- 
radial valence terms of increasing orders while re- adjusting p until convergence in the 
binding energy is reached, see Figure HI Starting with a converged wavefunction, different 
values of po are tested using correlated sampling to finally locate the energy minimum, as 
shown in Figure [5j The converged MN wavefunction contains all spin-singlet valence terms 
with K < 12 and rii ag < 5 in both Y and T Jacobi bases. For the MN-SO case, the spin- 
singlet and spin-triplet valence states are mixed, and after antisymmetrization with the 
core, many components become almost orthogonal which creates numerical noise. It was 
for this reason that a more advanced optimization technique — comparative optimization 
on two independent random walks — needed to be developed. The energy minimum with 
Po is first found for an auxiliary wavefunction as shown in Figure |5j and the MN-SO 
wavefunction is then tailored to the optimum value of po- Only valence terms lowering 
the energy the most are admitted to the MN-SO wavefunction at any optimization stage, 
or in other words n[ ag and all components of j va i in Eq. ( Ti~5"j) are treated as discrete 
variational parameters. It takes ninety carefully selected valence terms with K < 14 
and riiag < 5 to reach energy convergence in the MN-SO case (see Figure H]). Due to 
a more selective optimization process, the MN-SO wavefunction contains fewer valence 
terms than its MN counterpart. 

In both cases, the energy minimum in Figure is located at p ~ 0.45 fm. In either 
case, the optimization begins with a valence term having low hyper-momentum K and 
n iag = 0. Because of their diverging local kinetic energies near the origin, valence terms 
with ni a g 7^ are avoided until preliminary convergence with K has been reached. To 
avoid high partial waves in the valence part, both Y and T Jacobi configurations are mixed. 
The linear coefficients c in Eq. ([T]) are determined via energy matrix diagonalization. To 
report final results, values of observables obtained on several independent random walks 
are averaged for precision. Due to the Monte Carlo integrations, all results come with 
statistical errors. 

3.1. Three-body calculations 

For comparison purposes, we repeat the inert-core three-body calculations for 6 He pub- 
lished in [H2]- Here, we make use of the valence basis introduced in Section 12.21 The phe- 
nomenological 4 He-n interaction, vanishing for f- and higher-order partial waves, combines 
a spherically symmetric central Woods-Saxon and a spin-orbit Woods-Saxon-derivative 
parts with parameters taken from [ 42J. As for the interaction between valence neutrons, 
a realistic Gogny force is used [H3]. Although the core-n and n-n two-body interactions 
reproduce a-n and n-n phase shifts satisfactorily, 6 He bound by these interactions would 
miss more than half of its experimental three-body binding energy of -0.97 MeV. As in 
[ H2] , the problem of underbinding [ [2J H] commonly encountered in three-body models 
based on inert cores is cured with an effective three-body force presumably simulating the 
effects of the closed 3 H+ 3 H channel. 

Three-body calculations were performed using the computer code efadd [SI]. The 
Pauli principle was satisfied by writing the wavefunction in the T Jacobi basis (see 
Section |2T2"|) and by projecting the forbidden core-n states before diagonalization [[7]. 
Upon fitting the strength of the ambiguous three-body force, the converged three-body 
binding energy is -0.98 MeV for the wavefunction containing all valence terms with K < 40 
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Figure 4. Convergence of the three-body binding energy with the number of valence 
terms included in Eq. ([1]). In case MN, p Q is adjusted along the optimization route; in 
case MN-SO, the results are for fixed po — 0.45 fm. Error bars were not computed for all 
points, and even when present they may be smaller than actual symbols. 



and niag < 25. More three-body results are shown in Section H] alongside with those from 
our microscopic model developed in this work. 

In three-body models, the size of the core does not enter the actual calculations and 
some ad hoc assumption about the radius of the core is needed to estimate the size of the 
whole system. Then, the rms point proton (r 2 ) 1 ^ 2 and rms point nucleon matter (t" 2 ^) 1 / 2 
radii of the three-body system (mass number A) are related to those of the core (mass 
number A core ) through: 



K) 1 ' 2 = J( r2 P ^re)) + (r 



core-CMS 



>, (r 2 J 



2U/2 
ml 



A core (rl(core)) + (p 2 ) , (26) 



where r core _ C MS is the distance between the core's center of mass and the center of mass 
of the whole nucleus, and () denote expectation values. 



4. Results 

With the fully optimized wavef unctions, we calculate binding energies, rms radii, and 
density distributions, study core-valence antisymmetrization effects, and extract two- 
neutron overlap functions for 6 He. 

4.1. Energies and radii 

Binding energies and radii for 4 He and 6 He from our calculations are shown in Table Q] 
along with experimental values and results obtained in the three-body model described 

SVM 



in Section [371 
calculations 



19] representing microscopic cluster models, and ab-initio GFMC 



In Table [T], experimental rms point proton radii 
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Figure 5. Dependence of the three-body binding energy of converged 6 He on the non-linear 
parameter po- Both curves are constructed in correlated sampling on walks produced for 
reference values of p depicted by stars. In the MN-SO case, error bars could be computed 
for the reference po only 



from accurately measured charge radii (r 2 ) 1 ^ 2 by using the relationship [ |4"6 



(4) 



(r 2 c )-(R 2 p )-(Rl)^, 



(27) 



where (R 2 ) 1 ^ 2 = 0.895(18) fm [ H7] is the rms charge radius of the proton, (R 2 ,) = 
—0.120(5) fm 2 [ HS1 HH] is the mean-square charge radius of the neutron, and N and Z 
are the neutron and proton numbers, respectively. 

Based on arguments in Section [1] and Section I3.1[ there is a qualitative difference be- 
tween three-body and microscopic results in Table [TJ Strictly speaking, the three-body 
results should be taken with caution because the three-body binding energy was fitted 
using an auxiliary three-body force of an arbitrary strength, while for radii we arbitrarily 
assumed the size of the core in Eq. ( 12^|) to be equal to that of our MN-SO 4 He. The 
later assumption was made to ensure the best comparison between three-body and our 
results. Similar arbitrary assumptions are made for densities from the three-body model 
in Section I4T21 

In our model, a free 4 He in Table [T] is overbound and smaller relative to experimental 
data which sets a wrong scale for absolute binding of 6 He. We are, however, mostly inter- 
ested in three-body-like features of 6 He which should depend more on three-body rather 
than the absolute binding energy. As mentioned in Section 12. 3[ the mixture parameter in 
the Minnesota force was adjusted so that the three-body binding energy is about right. 
Upon this adjustment, there is no additional freedom in computing other observables. 

Most likely due to the smaller 4 He core, the MN and MN-SO proton radii of 6 He are 
smaller than they should be. On the other hand, matter radii are comparable with those 
deduced from experiments. To assess how strongly the radii of 6 He depend on the size of 
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Table 1 

Absolute (E) and three-body (E^ody) binding energies in [MeV], and rms point nucleon 
(proton p, neutron n, matter m) radii in [fm] of 4 He and 6 He. MN and MN-SO are 
results of this work; for other models, see text. Experimental proton radii were com- 
puted from Eq. (|2"7|) using charge radii from references cited in the table. Experimental 
rms neutron radii were computed from experimental values of proton and matter radii 

using (r^) = (l/A)[Z(r 2 ) + N(r 2 )]. The thickness of the neutron halo is defined as 
Ar=(r2) 1 /2_( r 2)i/2_ 
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the core, we turn to the three-body model. In a naive three-body picture based on Eq. (1261) , 
the larger size of 6 He in the three-body model (when compared to MN and MN-SO) is 
solely due to the valence neutrons living on average slightly farther from the core. If the 
radius of the core is increased to its experimental value 1.46 fm, the three-body radii of 
6 He become {r 2 m ) 1/2 = 2.51 fm, {r 2 p ) l l 2 = 1.90 fm, {r 2 n ) l l 2 = 2.77 fm, and Ar = 0.86 fm, 
and the experimental proton radius of 6 He would be seemingly reproduced. A similar 
shift towards larger radii could be expected for our results if larger cores were involved. 
Perhaps due to a stronger three-body binding, MN-SO 6 He is slightly smaller than its 
MN counterpart, but the thickness of the neutron halo remains about the same. 

Within SVM, 6 He has been studied in the past repeatedly, e.g. [[191 E5]- In Table [Q 
SVM results obtained in [ [19] for central and spin-orbit Minnesota and Coulomb interac- 
tions are quoted. In that reference, several different cluster compositions were considered 
to study break-up of the core in 6 He. For 4 He, we quote the results for the «2 model 
where the 4 He wavefunction is a superposition of three Os-harmonic oscillator Slater de- 
terminants with common oscillator parameters. Due to this simple picture, SVM 4 He is 
bound significantly less compared to MN and MN-SO cases. The SVM results for 6 He 
in Table [1] are those from model (b) in [ 19|. In that model, 6 He was described as a 
combination of 4 He+n-|-n and 3 H+ 3 H with tritons again built from simple Os-harmonic 
oscillators. The triton channel was introduced to overcome the insufficient three-body 
binding of 6 He. SVM radii and three-body binding energies of 6 He are comparable with 
ours, especially with the MN model. 
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Figure 6. Point proton and neutron densities in 6 He from models in Table [JJ (except 
MN). SVM and GFMC densities are from [[19] and [|9], respectively; for densities from 
the three-body model, see text. For comparison, the proton (=neutron) density of the 
MN-SO 4 He is shown. All proton (neutron) distributions are normalized to the number 
of protons (neutrons). Here, r is the distance from the center of mass. 



For the sake of completeness, we also show ab-initio GFMC results in Table [TJ These 
were obtained using realistic two-body AV18 and three-body IL2 interactions. The three- 
body binding energy for this case was computed by using S( 6 He) = —29.4(1) MeV | 
143] . GFMC results are shown to point out that, by using modern realistic potentials in 
microscopic calculations, absolute binding energies and proton radii of 4 He and 6 He can 
indeed be reproduced. However, as we argued in Section [1] and unless proven otherwise, 
questions may arise about how well ab-initio models treat asymptotic regions so important 
for Borromean halo nuclei. Also, two-neutron overlap functions are yet to be extracted 
from ab-initio models of two-neutron halo nuclei. 

4.2. Density distributions 

Point nucleon density distributions in 6 He for the more realistic MN-SO case are plot- 
ted in Figure [6] along with those obtained from other models listed in Table [TJ To com- 
pute densities from the three-body model, some assumption is needed about the internal 
structure of the core; quite arbitrarily, we computed these densities using an auxiliary 
microscopic wavefunction of the type of Eq. ([TJ constructed by combining the MN-SO 
4 He core and the three-body wavefunction serving as the valence part. Once again, this 
arbitrary construction reveals deficiencies in the few-body models. 

The densities in Figure [6] from different models are close to one another with small 
differences reflecting different radii and wavefunction compositions. All models reproduce 
the most pronounced property, the neutron halo with the neutron distribution extending 
far beyond that of protons. Depleted at short distances, the proton density of 6 He stretches 
farther out than that of 4 He. A partial explanation of this effect comes from the three- 
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body model: in 6 He, the a core does not sit at the center of mass of the entire system, 
and its motion relative to the center of mass spreads out the proton distribution. Due 
to the same effect, the neutron density in 6 He is also expected to be depleted at small 
distances relative to that of a free 4 He, as is also visible in Figure |6j 

4.3. Antisymmetrization effects 

An important drawback of three-body models is the approximate way in which the 
Pauli principle between the core and the valence particles is taken into account [|3]. The 
three-body procedure closest to the proper antisymmetrization is Feshbach projection 
where forbidden states are projected before diagonalization [[7J. Given the similarities 
between the three-body basis and our valence basis, the effects of Pauli blocking can 
now be assessed microscopically from our model. This can be done simply by including 
or not including the core- valence antisymmetrizer ¥ \ core ~ val i n Eq. ([T]). Because of the 
simpler optimization procedure involved, we performed this study for the MN 6 He only, 
but qualitatively the same outcome is also expected for the MN-SO case. See [ [27] for 
details. 

The valence channels with K(= l x = l 2 ) = suffer the most from the core-valence Pauli 
blocking. When j[ core - val [ s not included and the K = valence channels are present in 
the wavefunction, 6 He is three-body overbound by several tens of MeV. When all K = 
channels are removed, the nucleus becomes three-body unbound regardless of the inclusion 
of valence terms with higher hyper-momenta. When j\ core ~ val is active, converged 6 He 
including K = valence channels is three-body bound by about -0.9 MeV as shown in 
Table [TJ Upon removal of all K = channels from the converged wavefunction, the three- 
body binding reduces to about -0.75 MeV. 

It is evident that the proper antisymmetrization is crucial for the structure of 6 He. To 
produce a meaningful 6 He it is not sufficient to simply neglect the most Pauli-blocked 
K = valence channels. Rather, all contributing valence channels ought to be included 
in the model space and carefully antisymmetrized. Given this conclusion, we strongly 
advocate using the best available Pauli blocking techniques in three-body models of halo 
nuclei. 

4.4. Overlap functions 

All models mentioned in Table [U although different in their nature and predictive 
power, are in fair agreement on the most commonly computed properties of 6 He. To 
appreciate the amount of microscopic details embedded in different models, it would be 
better to compare these models at the level of wavefunctions rather than highly integrated 
observables. In SVM, only s-wave overlap functions for 6 He have been computed, but not 
expanded in hyper-spherical coordinates [US]; in GFMC, these functions are yet to be 
computed. 

Here, we compare overlap functions (Section 12. 4p computed for the MN-SO 6 He to 
three-body wavefunctions (Section l3.ip . In both cases, the 4 He core is in its ground state. 
The three-body wavefunction is normalized to unity. For a meaningful interpretation of 
overlap functions, MN-SO wavefunctions of 4 He and 6 He in Eq. (T2"0"j) need to be normalized 
to unity. For 4 He, the normalization is known analytically from SVM; the norm of the 
6 He wavefunction was determined numerically with accuracy of 0.3% or better by using 
an auxiliary sampling function [ [27J . 
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Table 2 

Spectroscopic factors of the five dominant overlap channels in 6 He. All channels are in 

the T Jacobi basis. Numbers in parentheses are relative errors. 
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Ordered by spectroscopic factors, the five strongest overlap channels in the MN-SO 
6 He are listed in Table [2j These are the only channels that could be resolved, all other 
potential channels have spectroscopic factors too small and as such are buried in numerical 
noise. The table also contains three-body results; these were obtained by summing up all 
hyper-radial components in a given three-body channel. Not only the dominant channels 
are the same in the two models, but also the order of their spectroscopic strength is 
preserved. In the three-body model, these five channels account for more than 98% of 
the wavefunction. Therefore, we expect that these channels should also grasp most of the 
4 He+n+n decomposition of the MN-SO 6 He ground state. 

The hyper-radial dependence of overlap channels from Table [2] is shown in Figure [7J 
and Figure El We show these functions as u(p) because of their simpler asymptotic fall- 
off presented in Eq. ( 12"B"|) . The three-body and our model agree on the overall shape of 
their hyper-radial distributions, and the number of nodes in those channels. On the other 
hand, our model appears to narrow the hyper-radial distributions and shift their peaks 
towards smaller distances, and to put more weight on smaller hyper-radii. This difference 
may be related to smaller radii of the MN-SO 6 He in Tabled] when compared to those of 
the three-body model. 

The most striking difference between our overlap functions and three-body wavefunc- 
tions is the normalization of their components quantified in Table [2] and visible in Figure [7J 
and Figure [H The MN-SO spectroscopic factor for the strongest K = 2 s-waves channel is 
greater than one due to recoil effects as the 4 He core in Eq. ffTSl) does not sit at the center 
of mass of 6 He. These recoil effects are present in the three-body model too, but there the 
wavefunction as a whole is normalized to unity and so all spectroscopic factors are smaller 
than one. This deficiency of three-body models has been pointed out in [ [5] where an up- 
per limit 25/16=1.5625 was derived for an additional renormalization factor to multiply 
the three-body spectroscopic factor as a mock-up for missing microscopic information. 
Indeed, our microscopic model consistently predicts spectroscopic factors greater by at 
least 30% than those obtained in the three-body model, but the increase varies between 
channels in a non-trivial manner as can be seen in Table [2j This observation suggests 
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Figure 7. Hyper-radial dependence of overlap MN-SO functions and three-body wave- 
functions for 6 He for those channels presented in Table [2j The legend is the same in all 
panels. 



that to account for microscopic effects, it may not be sufficient to simply renormalize the 
entire three-body wavefunction by a common factor (such as a/25/16 suggested in [[5]). 
This conclusion is important, for example, for two-neutron transfer reaction theories for 
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6 He in which three-body wavefunctions are traditionally used as structure input. 

Since two-neutron transfer cross sections depend on spectroscopic factors, using micro- 
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scopically derived overlap functions instead of three-body wavefunctions would certainly 
have implications for reaction observables. For example, if the transfer reaction is mostly 
sensitive to the dominant component (s-waves) and peripheral, one might not see much 
change in the normalization, since the asymptotic parts of the overlaps do not signifi- 
cantly change. In the other hand, if the reaction is sensitive to the whole volume of the 
nucleus, one can expect a renormalization of the cross section consistent with the spectro- 
scopic factors. Often, two-neutron transfer reactions are complicated by interference of 
various mechanisms. Preliminary two-nucleon transfer calculations for 6 He(p,t) 4 He at 
Ei a b = 25 MeV involving the 6 He microscopic overlap functions here presented have been 
performed assuming a 1-step reaction. However, any meaningful comparison with the 
data requires an extended study of the reaction mechanism which falls beyond the scope 
of this work. 

At hyper-radii beyond about 12 fm, the MN-SO overlap functions in Figure [7J and 
Figure M become unreliable due to statistical fluctuations. Very large hyper-radii would 
place valence neutrons into regions very distant from the core, and because the Monte 
Carlo sampling probability is proportional to the wavefunction squared, such extreme 
spatial configurations are very unlikely to be visited by a walker during a random walk. 
Moreover, statistical samples in such distant regions may be highly correlated. For exam- 
ple, in extreme configurations of 6 He, a hyper-radius of 12 fm corresponds to a di-neutron 
at a distance of 10.4 fm from the center of the core, or to two neutrons on opposite sides 
of the core, 17 fm apart. 

In asymptotic regions, where the two valence neutrons are distant from the core, the 
core-valence antisymmetrization effects in Eq. ([T]) disappear. Then, both the three-body 
wavefunction and the overlap functions should fall off exponentially as shown in Eq. (^ 
and Eq. (1231) . For E^ody ~ — 1 MeV, we get k ~ 0.22 fm -1 . At larger hyper-radii in 
Figure El the three-body wavefunctions have the right slope still influenced by the three- 
body centrifugal barrier. Interestingly, we found in Section [3] the value po — 0.45 fm 
to be optimal for the valence part of the fully antisymmetrized wavefunction of 6 He. By 
comparing Eq. (fT3"j) with Eq. (12"31 . the two decay parameters can be related as k = l/(2p ). 
Then, po — 0.45 fm would correspond to k = 1.11 fm -1 , or Embody ~ —25 MeV. In 
other words, all individual valence terms in the MN-SO wavefunction decay much faster 
than the expected asymptotic form of overlap functions. Nevertheless, the strongest 
MN-SO overlap functions follow approximately the right long-range trend within the 
computationally safe region. Albeit the efforts in tailoring the basis to the problem, 
our results still show lingering differences in the asymptotics which could be improved 
by inclusion of Laguerre polynomials of even higher orders and/or by using a value of 
Po > 0.45 fm in the MN-SO wavefunction. 

Even though not discussed in detail, we also computed overlap functions for the MN case 
of 6 He. Due to the absence of the spin-orbit force, valence spin-triplet states are missing 
in the MN wavefunction and consequently among overlap functions. Spectroscopic factors 
of spin-singlet channels were about the same as those for the MN-SO model in Table [2j 

4.5. Clusterization 

Although 6 He has been the object of many experimental studies (e.g. [[561 [57]), there 
is still an ongoing debate on the neutron clusterization in this nucleus, namely on the 
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Figure 9. Correlation density plot for the ground state of 6 He using MN-SO. Two domi- 
nant clustering patterns are shown schematically. 



dominance of the cigar- or the di-neutron-like configuration. One way of gaining insight 
into this is by calculating the probability of finding the valence neutrons at definite places 
within the three-body- like decomposition defined in Eq. (T2~Tj) : 
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where r n _ n and r core _ nn are the distance between valence neutrons and the distance be- 
tween centers of masses of the core and the valence neutron pair, respectively. The 
integration in Eq. ( 128]) is carried out over all (unprimed) degrees of freedom of all nu- 
cleons. The probability plot for the MN-SO 6 He is presented in Figure |9j The figure 
indeed exhibits two peaks, as it should for a system dominated by K = 2 s— waves in the 
T— basis: a di-neutron-like peak positioned at about r n _ n = 1.93 fm and r core _ nn = 2.63 fm 
(p = 3.33 fm) with the two neutrons close together and far from 4 He, and a cigar-like 
peak at r n _ n = 3.82 fm and r core _ nn = 1.03 fm (p = 2.95 fm) with the two neutrons po- 
sitioned on opposite sides of the core. Qualitatively the same clustering picture has been 
predicted for example within three-body models [[2] and SVM [ES]. We looked into the 
clusterization probability distribution produced within the three-body model and found 
the peaks to be shifted to slightly larger radii (the di-neutron-like peak at r n _ n = 2.08 fm 
and r COT . e _ nn = 2.88 fm and the cigar-like peak at r n _ n = 4.18 fm and r core _ nn = 1.08 fm) 
due to the three-body wavefunctions peaking at larger hyper-radii in Figure [71 Integrals 
under the peaks in the microscopic model are larger than those in the three-body model, 
consistent with the spectroscopic factors presented in Table HJ the integral under the 
di-neutron peak for the microscopic calculation is 0.85, compared to 0.60 in the three- 
body model and the integral under the cigar-like peak for the microscopic calculation is 
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0.54, compared to 0.38 in the three-body model. In both models, it appears that there 
is coexistence of the two clustering patterns, with 60% di-neutron and 40% cigar-like 
configurations. 



5. Summary and Outlook 

In this work, a microscopic structure model for light two-neutron halo nuclei is devel- 
oped. Our goal was to combine advantages of few-body and microscopic methods to create 
a model capable to deal simultaneously with short- and characteristic long-range halo ef- 
fects exhibited by these nuclei. To succeed, we combine the stochastic variational method 
and the hyper-spherical harmonic method into a fully antisymmetrized many-body ap- 
proach. From the computational point of view, our model is similar to the variational 
Monte Carlo method, the novel feature being the form of the basis which incorporates 
the few-body features of two-neutron halo nuclei, in particular the correct behavior when 
the two halo neutrons are far from the core. 

The model is applied to the ground state of 6 He bound by an effective Minnesota 
nucleon-nucleon interaction including the spin-orbit force. When comparing three-body 
binding energies, radii, and densities, our results are comparable with those from other 
microscopic models using the same interaction. The overall binding energy is not repro- 
duced as that would require more realistic forces, but this is not essential to produce the 
characteristic halo features determined by the binding relative to the three-body 4 He+n+n 
threshold. The halo nature of the nucleus can be seen from its extended neutron density 
resulting in a large difference between the rms matter and proton radii. We advocate 
that the standard highly integrated observables, such as binding energies and radii, are 
not sensitive to details of the halo part of configuration space. To recognize these details, 
one should compare structure models at the level of wavefunctions or look at observables 
more sensitive to long-range correlations, such as reaction observables. 

We properly calculate, to our knowledge for the first time, microscopic two-neutron 
overlap functions for 6 He and compare them with those from the state-of-the-art three- 
body model. These overlap functions provide a crucial input to reaction calculations 
involving 6 He, in particular to two-neutron transfer reaction models. As our basis was 
tailored specifically to capture the asymptotic behavior of the wavefunction, it is not sur- 
prising that the three-body wavefunction and the overlap functions are indeed similar in 
this region of coordinate space. In the range of the nuclear interaction, the three-body 
wavefunction seems to reproduce the properties of the many-body overlap functions only 
qualitatively. At a quantitative level, however, there are significant differences between 
overlap functions and three-body wavefunctions. Outstanding amongst these, the spectro- 
scopic factor for the dominant microscopic overlap channel is larger by about 40% than 
its three-body counter-part. This discrepancy reveals the deficiency of three-body mod- 
els, namely their inert-core approximation. Also, our study shows that the underbinding 
problem in three-body models comes partly from a poor treatment of antisymmetrization. 

It is often thought that few-body models can be corrected for many-body effects by a 
simple renormalization of the wavefunction. However, we demonstrate that, in general, 
such a simple fix may not be sufficient because for the most important three-body-like 
components in 6 He, the ratio between microscopic and three-body spectroscopic factors 
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varies between 1.3 and 1.7 in a nontrivial manner. 

In this work we have shown that it is possible to obtain the correct long-range properties 
of halo nuclei within a fully microscopic approach. This opens up the opportunity to 
many applications. Studies of two-neutron transfer reactions involving 6 He using our 
two-neutron overlap functions are underway. Given that in our model all matrix elements 
are calculated numerically, an extension of the method towards more realistic nuclear 
interactions should be straightforward. It would also be exciting to apply this method to 
heavier systems such as n Li which would most likely entail more advanced optimization 
algorithms. On top of that, further code parallelization would be necessary to tackle the 
factorial growth in dimensionality of the problem. 
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